Scaling of the specific heat in superfiuid films 
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We study the specific heat of the x — y model on lattices L x L x H with L ^> H (i.e. on lattices 
representing a film geometry) using the Cluster Monte-Carlo method. In the H -direction we apply 
Dirichlet boundary conditions so that the order parameter in the top and bottom layers is zero. We 
find that our results for the specific heat of various thickness size H collapse on the same universal 
scaling function. The extracted scaling function of the specific heat is in good agreement with the 
experimentally determined universal scaling function using no free parameters. 
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The theory of second order phase transitions is based 
on the assumption that at temperatures close to the crit- 
ical temperature T c there is only one dominating length 
scale associated with the critical behavior of the system, 
the correlation length. Because the correlation length 
diverges as the critical temperature is approached the 
microscopic details of the system become irrelevant for 
the exponents describing the most singular dependence 
of the thermodynamic functions on the reduced temper- 
ature t = T/T c — 1. This intuitive picture has its founda- 
tion in the renormalization group treatment of second or- 
der phase transitions. Within the renormalization group 
treatment it becomes evident that the critical behavior 
can be divided into different universality classes which 
are characterized by a set of critical exponents. 

If the system is confined in a finite geometry (e.g. a 
cubic or film geometry) the finite-size scaling theory [1] 
is thought to describe well the behavior of the system 
at temperatures near T c . The intuitive idea behind the 
finite-size scaling theory is that finite-size effects are ob- 
served when the bulk correlation length £ becomes of the 
order of the finite system size, i.e. for our case here the 
film thickness H . For a physical quantity O this state- 
ment can be expressed as follows [2]: 
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The dimensionless ratio on the left-hand-side of the above 
equation is a universal function f(x) of the dimensionless 
ratio x = H/£, i.e. in our present work the film thickness 
measured with respect to the correlation length. The 
function / depends on the boundary conditions and the 
geometry of the system. 

Liquid helium 4 He has been a very good ground both 
for testing finite-size scaling theory and measuring the 
critical exponents that go along with a second order 
phase transition in the case of a complex order param- 
eter. Measurements of the superfiuid density [3] and 
specific heat [4] on helium films, however, fail to ver- 
ify the finite-size scaling theory. Field theoretical cal- 
culations for the standard Landau-Ginzburg free energy 
functional in different geometries with Dirichlet bound- 
ary conditions have been carried out [5-8]. New specific 



heat measurements [9] and also a reanalysis [10] of the 
old specific heat data [4] show good agreement between 
the results of the calculations reported in Refs. [5-7] and 
those data. Furthermore, new experiments on liquid 4 He 
under microgravity conditions are planned [11] to exam- 
ine the finite-size scaling properties of the specific heat. 
The above mentioned field theoretical calculations are 
based on methods such as a loop expansion which system- 
atically includes corrections due to fluctuations around a 
mean-field solution. These methods are approximate in 
nature and assume that the basic physics of the mean 
field treatment is qualitatively correct. In addition these 
methods neglect the effect of vortices. Thus, accurate nu- 
merical calculations of the specific heat of helium films 
and its scaling with film thickness are welcome. 

In our previous work, we used the x — y model with 
periodic boundary conditions [12,13] in the direction of 
the film thickness H to compute the superfiuid density 
and the specific heat of thin helium films. We demon- 
strated scaling with the expected values for the critical 
exponents of the superfiuid density and the specific heat 
with respect to the film thickness, thus confirming the 
validity of the finite-size scaling theory. However, the 
obtained universal function for the specific heat does not 
match the experimentally determined universal function 
of Ref. [9], indicating that periodic boundary conditions 
are only a poor approximation of the correct physical 
boundary conditions. 

In this paper we study the effect of Dirichlet bound- 
ary conditions on the finite-size scaling behavior of the 
specific heat of 4 He in a film geometry L x L x H with 
L H . These boundary conditions are believed to sim- 
ulate the boundary conditions of the real helium films 
more accurately [6,7]. With these boundary conditions, 
the superfiuid order parameter vanishes in the top and 
the bottom of the film. The specific heat of helium films 
has been measured in films with very large planar di- 
mensions. The large planar dimension (L) of the films is 
taken into account in our calculations by assuming peri- 
odic boundary conditions in the planar ^-directions and 
L H so that the finite-size effects on our calculated 
specific heat due to L are not important. We find that 
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our results for the specific heat obey the finite-size scal- 
ing theory using the expected values of the critical ex- 
ponent v. More importantly, our calculated universal 
scaling function agrees well with that obtained by the re- 
cent measuments of the specific heat of helium films [9] 
with no adjustable parameters. 

We use the x — y model to describe the fluctuations of 
the order parameter in superfluid 4 He near the A-critical 
point (cf. e.g. Ref. [14]). In the pseudospin notation, the 
x — y model takes the following form: 
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where the summation is over all nearest neighbors, s = 
(cos 9, sin 9), and J sets the energy scale. The angle 9 
corresponds to the phase of the superfluid order param- 
eter t/>(r). This order parameter is the average value of 
the helium atom creation operator which is defined in a 
volume whose linear extensions are much larger than the 
interparticle spacing and much smaller than the corre- 
lation length. This condition can be realized only very 
near the critical temperature. 

In our calculations Dirichlet boundary conditions can 
be imposed by coupling the top and bottom layers of 
the film to two static staggered pseudospin configurations 
playing the role of the "substrate" layers so that the pseu- 
domagnetization (the superfluid order parameter) in the 
top and bottom substrates is exactly zero, i.e. 

s(x,y,z) = (-iy+ys(l,l,z), z = 0,H + l, (3) 

where x and y label the integer coordinates of the lat- 
tice sites in the two "substrate" planes perpendicular to 
the if-direction one below the film, i.e. z = 0, and one 
above the film, i.e. z = H + 1. In the ^-directions we 
applied periodic boudary conditions. The crucial differ- 
ence between the boundary conditions used in this work 
and periodic boundary conditions in all directions [12,13] 
is that the superfluid density develops a profile in the 
if-direction, whereas it is completely homogeneous for 
periodic boundary conditions. 

We computed the specific heat on L 2 x H size lattices, 
where L = 20, 40, 60, 80, 100 and H = 12, 16, 20, 24. The 
specific heat c is obtained by 



c=t( { H 2 )-{Hf) 



(4) 



where [3 = 1/(^7) and N is the number of pseudospins 
contributing to the specific heat. The multi-dimensional 
integrals in the expression for the average is computed by 
means of the Monte-Carlo (MC) method using Wolff's 
cluster algorithm [19]. We carried out of the order of 
20, 000 thermalization steps and of the order of 750, 000 
measurements. The calculations were performed on a 
heterogeneous environment of computers including Sun, 
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FIG. 1. The specific heat as a function of T for L 2 x 12 
lattices. In the 3D x - y model T x /J = 2.2017 [17]. 



IBM RS/6000 and DEC alpha AXP workstations and a 
Cray-YMP. 

Since we do not possess an easy procedure to extrapo- 
late to the L — ► oo limit for the values of the specific heat 
computed on finite lattices L x L x H , we shall use our 
largest value of L = 100 to represent the infinite L limit. 
This can be justified because the specific heat appears 
independent of L for L > 60 (cf. Fig.l). Furthermore, 
for given thickness H < 24 used in this work, we do not 
expect the maximum of the specific heat to grow signif- 
icantly with increasing values of L > 100 because for 
temperatures in the range T 2D (H) < T < T\ (where 
T 2D (H) is the Kosterlitz-Thouless critical temperature 
for films of thickness H) the behavior of the specific heat 
can be described by Kosterlitz-Thouless theory which 
leads to a finite value of this maximum. In order to il- 
lustrate this argument further, we show in Fig. 2 the size 
dependence of the specific heat c(T, L) computed on pure 
two-dimensional lattices L x L (up to L = 400) with pe- 
riodic boundary conditions. The Independence of the 
specific heat can be neglected for values of L > 80. 

Now, we would like to check the finite-size scaling hy- 
pothesis for the specific heat c(t, H) for large L (L H). 
We consider the finite-size scaling expression for the spe- 
cific heat c given by [5,6]: 



c(t,H) = c(t ,oo) + H a ' v hitH 1 '"). 



(5) 



The function fi(x) is universal and v = 0.6705 as has 
been extracted from recent accurate experiments [16]. 
The hyperscaling relation a = 2 — 3^ yields ajv = 
—0.0172. At the reduced temperature to the correlation 
length £(t) = \t \~ v becomes equal to the film thickness 
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FIG. 2. The specific heat for pure two-dimensional lattices 
L x L with periodic boundary conditions. 

H, i.e. t = (Co/H) 1 ^ with ^+ = 0.498 [18]. We have 

c(i ,oo) = c(0,oo) + 2+*- a , (6) 

where we use the bulk values c(0,oo) = 30, = —30 
obtained by studying the finite-size scaling of the specific 
heat of cubes [13]. 

Thus, we can make a direct comparison of our universal 
function fi(x) to the experimentally determined scaling 
function fi(x) given in Refs. [9]. In order to do that we 
express all lattice units in physical units leading to the 
conversion formula: 

M X )\phy, = M X )Lttice ( 7 ) 

where V m is the molar volume of 4 He at saturated va- 
por pressure at T\. The unit of length a (i.e., the lat- 
tice spacing a in the x — y model) can be determined 
in such a way that the x — y model describes the fluc- 
tuations of the order parameter in liquid helium. This 
means that the calculated bulk helicity modulus T from 
the x — y model (which is related to the superfluid den- 
sity p s by p s = (m/h) 2 T where m is the bare 4 He atom 
mass) and the one measured experimentally match near 
the critical point. Namely, /?T has dimensions of inverse 
length and near T\ can be expressed in units of a as 
f3T(t — ► 0) = Yot v I a which is matched with the experi- 
mental temperature dependence to obtain the value of a 
in Angstroms. We find a = 2.95A [13] and thus 

/i(*)U. = 15 - 02 /i(*)L«.-ce Joule/(°Kmole). (8) 

Therefore, the function fi(x) in physical units is com- 
pletely determined. 



FIG. 3. The universal function fi(x) for films with Dirich- 
let boundary conditions is compared to the results of our ear- 
lier calculation using periodic boundary conditions [13] (open 
squares). The solid lines are the results to our fit to the func- 
tion given by Eq. (10). The dashed line represents the results 
of the loop expansion calculation [5,6]. 

In Fig. 3 we give the results of our present Monte Carlo 
calculation with Dirichlet boundary conditions along the 
film thickness direction and we compare them to the re- 
sults of our earlier Monte Carlo calculation [13] using 
periodic boundary conditions. Each set of our Monte 
Carlo data collapse onto a universal curve which defines 
the function fi(x) for periodic and a different function 
fi(x) for Dirichlet boundary conditions. As a compari- 
son we also plot the result for fi(x) of the field theoreti- 
cal calculation for the Landau-Ginzburg functional in a 
film geometry with Dirichlet boundary conditions [5,6] 
(dashed line). 

The universal function fi(x) satisfies the following lim- 
its 

fl ( x ^ ±00 )^ c± \ x \-», ( 9 ) 

For sufficiently large values of \x\ the universal func- 
tion fi(x) for the film geometry should agree with 
the universal function obtained for the cubic geome- 
try (bulk). In our earlier work [13] we found that 
for cubes c + = — 449. 6 Joules /(° K mole) and c_ = 
—430.8 J oules / (" Kmole) for a = —0.0115. Using these 
asymptotic constraints on the form of the function fi(x) 
we were able to fit our results for fi(x) for intermediate 
and large values of |a? | with the formula 

f 1 (x) = c° ± + c ± \x + x° ± \- a . (10) 

where the + and — refer to the form which should be 
used for values of x above and below the peak of the 
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function fi(x) respectively. By fitting our results to the 
above function using the above mentioned values of c± 
we find that c° = 452.3 ± 0.3 Jou/es/(° Kmole) and c°_ = 



450.8 ± 0.2Joules/(°Kmole) and x° 



4.8 ± OAA and 



x_ = 3.4 ± 0.5A. Notice that the asymptotic form of 
Eq. (10) and the limits given by Eq. (9) agree. The 
asymptotic forms (10) describe our data accurately for 
-40 < x < -16 and -2 < x < 10. 
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FIG. 4. The universal function fi(x) for films with Dirich- 
let boundary conditions is compared to the experimental re- 
sults [9]. 

In Fig. 4 we show fi(x) obtained from our Monte-Carlo 
calculation (solid circles) using Dirichlet boundary condi- 
tions and the experimental data points taken from Ref. 
[9]. The agreement between theory and experiment is 
quite satisfactory. 

The range of the variable x = tH~ , in which our cal- 
culation of the scaling function fi(x) is valid, is deter- 
mined by the requirement £ a (the lattice spacing) for 
critical (continuum limit) behavior, which implies that 
t <C ([£o])" (Ko] is the prefactor of the correlation length 
= a[^o]t~" near the lambda point in lattice spacing 
units) and therefore \x\ <C ([£o]H)i . Taking the numer- 
ical values for [£o] ~ 1.2 and H ~ 20 used in our calcu- 
lations we find that we can expect scaling for \x\ <C 100. 
Indeed we find that for the values of H used in our cal- 
culations we obtain scaling for \x\ < 40. For \x\ > 40 
we find deviations from scaling. If we wish to obtain 
the scaling function for even larger values of \x\, we need 
to calculate the specific heat on even larger thickness 
helium films. However, the calculation of the specific 
heat for several temperature values on lattices of sizes 
100 x 100 x H with a given value of H, e.g. H = 24, us- 
ing cluster Monte Carlo, takes approximately 200 hours 



of Cray-YMP CPU time; thus, calculations for values of 
H much larger than the currently studied are beyond 
realistic computational time scales. 

The results of our calculations can be used in a 
straightforward way to deduce the specific heat values for 
any size helium films. Such a simple calculation might be 
useful for the preparation of the ground-based work for 
the confined helium experiment (CHeX) which is planned 
to be conducted in 1997 in space under microgravity con- 
ditions [11] where the rounding of the transition due to 
gravity can be avoided. 

This work was supported by the National Aeronautics 
and Space Administration under grant no. NAGW-3326. 



[1] M. E. Fisher and M. N. Barber, Phys. Rev. Lett. 28 1516 
(1972); M. E. Fisher, Rev. Mod. Phys. 46 597 (1974); V. 
Privman, Finite Size Scaling and Numerical Simulation 
of Statistical systems, Singapore: World Scientific 1990. 

[2] E. Brezin, J. Physique 43 15 (1982). 

[3] I. Rhee, F. M. Gasparini, D. J. Bishop, Phys. Rev. Lett. 

63, 410 (1989). 
[4] T. Chen and F. M. Gasparini, Phys. Rev. Lett. 40, 331 

(1978); F. M. Gasparini, T. Chen and B. Bhattacharrya, 

Phys. Rev. B 23, 5797 (1981). 
[5] R. Schmolke, A. Wacker, V. Dohm, and D. Frank, Phys- 

ica B165 & 166 575 (1990). 
[6] V. Dohm, Physica Scripta T49 46 (1993). 
[7] P. Sutter and V. Dohm, Physica B194-196 613 (1994); 

W. Huhn and V. Dohm, Phys. Rev. Lett. 61 1368 (1988). 
[8] M. Krech and S. Dietrich, Phys. Rev. A46 1886 (1992), 

1922 (1992). 

[9] J. A. Nissen, T. C. P. Chui, and J. A. Lipa, J. Low Temp. 
Phys. 92 353 (1993), Physica B194-196 615 (1994). 
[10] A. Wacker and V. Dohm, Physica B194-196 611 (1994). 
[11] J. A. Lipa, private communications. 

[12] N. Schultka and E. Manousakis, to be published in Phys. 
Rev. B; preprint cond-mat/9406014 (1994). 

[13] N. Schultka and E. Manousakis, preprint cond- 
mat/9502062 (1995). 

[14] H. Kleinert, Gauge Fields in Condensed Matter, Singa- 
pore: World Scientific 1989. 

[15] S. Dasgupta, D. Stauffer, and V. Dohm, Physica A213 
368 (1995). 

[16] L. S. Goldner and G. Ahlers, Phys. Rev. B45 13129 

(1992) . 

[17] W. Janke, Phys. Lett. A148 306 (1992). 

[18] A. P. Gottlob and M. Hasenbusch, Physica A201 593 

(1993) . 

[19] U. Wolff, Phys. Rev. Lett. 62 361 (1989). 



4 



